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ABSTRACT 

It is a well-established result that the propagation of linear transverse waves along a thin but isothermal 
magnetic flux tube is affected by the existence of the global cutoff frequency, which separates the propagating 
and non-propagating waves. In this paper, the wave propagation along a thin and non-isothermal flux tube 
is considered and a local cutoff frequency is derived. The effects of different temperature profiles on this 
local cutoff frequency are studied by considering different power-law temperature distributions as well as the 
semi-empirical VAL C model of the solar atmosphere. The obtained results show that the conditions for wave 
propagation strongly depend on the temperature gradients. Moreover, the local cutoff frequency calculated for 
the VAL C model gives constraints on the range of wave frequencies that are propagating in different parts 
of the solar atmosphere. These theoretically predicted constraints are compared to observational data and are 
used to discuss the role played by transverse tube waves in the atmospheric heating and dynamics, and in the 
excitation of solar atmospheric oscillations. 

Subject headings: magnetohydrodynamics (MHD) - Sun: atmosphere - waves 



1. INTRODUCTION 

Magnetic flux tubes existing in the photosphere and lower 
chromosphere of the Sun are considered to be narrow bun- 
dles of strong magnetic field lines that rapidly expand with 
height in the solar atmosphere (e.g., Solanki 1993). The 
fundamental modes of linear oscillations of these flux tubes 
are typically identified with longitudinal, transverse and tor- 
sional tube waves (e.g., Defouw 1976; Roberts & Webb 1978; 
Roberts 1979, 1981; Spruit 1981, 1982; Priest 1982; Hollweg 
1985; Roberts 1991; Roberts & Ulmschneider 1997), with the 
two latter waves being Alfven-like waves. 

Observational evidence for the existence of Alfven-like 
waves in different regions of the solar atmosphere was given 
by high-resolution observations performed by the Solar Opti- 
cal Telescope (SOT) and the X-Ray Telescope (XRT) onboard 
the Hinode Solar Observatory. According to De Pontieu et al. 
(2007) and Cirtain et al. (2007), signatures of Alfven waves 
were observed by the SOT and XRT instruments, respectively. 
Moreover, Alfven waves were also reported by Tomczyk et al. 
(2007), who used the Coronal Multi-Channel Polarimeter of 
the National Solar Observatory. Interpretations of these ob- 
servations were given by Van Doorsselaere et al. (2008) and 
Antolin et al. (2009), who concluded that the reported obser- 
vational results describe kink waves. As of today, the most 
convincing observational evidence for the existence of Alfven 
waves in the quiescent solar atmosphere was given by Mcin- 
tosh et al. (2011), who reported indirect evidence for such 
waves found in observations by Solar Dynamic Observatory 
(SDO). 

Moreover, Fujimura & Tsuneta (2009) used SOT observa- 
tions to study fluctuations in pores and intergranular magnetic 
structures, and concluded that such fluctuations could be ex- 
plained by the existence of longitudinal (sausage) and trans- 
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verse (kink) waves propagating along magnetic flux tubes em- 
bedded in the solar photosphere. They found oscillation peri- 
ods of 3-6 min for the pores and 4-9 min for the intergranular 
magnetic elements. More recently, Okamoto and De Pontieu 

(2011) used data from SOT and presented observational ev- 
idence for the existence of high-frequency transverse waves 
propagating along spiclues in the solar atmosphere. Similar 
high-frequency tranverse waves were reported by Yurchyshy n 
et al. (2012), who observed type II spicules using the New 
Solar Telescope. 

There is also observational evidence for the existence of 
torsional Alfven waves in the solar atmosphere as reported by 
Jess et al. (2009), who interpreted data obtained with high spa- 
tial resolution by the Swedish Solar Telescope (SST). Alfven- 
like motions were also observed by Bonet et al. (2008), who 
found vortex motions of G band bright points around down- 
flow zones in the solar photosphere, and by Wedemeyer- 
Bohm & Rouppe van der Voort (2009), who used data from 
the SST to demonstrate that more disorganized relative mo- 
tions of photospheric bright points can also induce swirl-like 
motions in the solar chromosphere. More recently, observa- 
tional evidence for Alfven-like waves was found by Bonet et 
al. (2010), who used the ballon-borne Sunrise Telescope, and 
by De Pontieu et al. (2012), who used the SST to find tor- 
sional motions in spicules. Moreover, magnetic swirls in the 
solar atmosphere were reported by Wedemeyer-Bohm et al. 

(2012) . 

Among the above described observational results, the ob- 
servations perfomed by Fujimura & Tsuneta (2009), Okamoto 
and De Pontieu (2011), and Yurchyshyn et al. (2012) are the 
most relevant because the authors directly refer to transverse 
tube waves, which are the main topic of this paper. Our 
goal here is to determine the propagation conditions for these 
waves by calculating the so-called cutoff frequencies. In gen- 
eral, there are global cutoff frequencies, which are the same 
along the entire length of the tube, and local cutoff frequen- 



2 



Routh, Musielak, & Hammer 



cies that are height dependent. If, under certain simplifying 
approximations, the cutoff frequency is constant along the en- 
tire length of the tube, it is called a global cutoff. Otherwise, 
if it is height dependent, it is a local cutoff. In Section 6, we 
compare these cutoff frequencies to the observational results. 

The global cutoff frequency represents the natural fre- 
quency of linear oscillations of the magnetic flux tubes, and its 
value restricts the wave propagation to only those frequencies 
that are higher than this cutoff. For isothermal and thin flux 
tubes, the global cutoff frequencies for longitudinal and trans- 
verse tube waves were first determined by Defouw (1976) and 
Spruit (1981), respectively. These global cutoff frequencies 
are ratios of the characteristic wave speeds to the pressure 
(density) scale height. Since the scale height and the wave 
speeds are constant along isothermal and exponentially ex- 
panding magnetic flux tubes, the resulting cutoffs are the same 
along the entire length of the tubes. The fact that there is no 
global cutoff frequency for torsional waves propagating along 
isothermal and thin magnetic flux tubes was demonstrated by 
Musielak, Routh & Hammer (2007). 

A method to derive the global cutoff frequency for longi- 
tudinal tube waves was introduced by Rae & Roberts (1982) 
and Musielak et al. (1987, 1995), who demonstrated that the 
wave equation for these waves can be transformed into its 
standard form (also referred to as the Klein-Gordon equation), 
which directly displays the global cutoff frequency. A sim- 
ilar method was used by Musielak & Ulmschneider (2001, 
2003) to determine the global cutoff frequency for transverse 
tube waves. These global cutoff frequencies are important 
in studies of atmospheric oscillations (e.g., Roberts 1991; 
Hasan & Kalkofen 1999; Hasan 2003; Musielak & Ulm- 
schneider 2003; Hasan 2008) and chromospheric and coronal 
heating (e.g., Narain & Ulmschneider 1996; Ulmschneider & 
Musielak 2003). 

The existence of global cutoff frequencies is restricted 
to thin and isothermal magnetic flux tubes. However, in 
more general cases when the flux tubes are either wide and 
isothermal, or thin and non-isothermal, or wide and non- 
isothermal, the resulting cutoff frequencies depend on at- 
mospheric height, which means that they are local quanti- 
ties (e.g., Routh, Musielak & Hammer 2007). A method to 
determine such local cutoff frequencies was developed by 
Musielak, Musielak & Mobashi (2006). This method re- 
quires that wave equations are cast in their standard forms 
(i.e., without terms with the first-order derivatives) and that the 
oscillation theorem (e.g., Kahn 1990) is used to obtain the cut- 
off frequencies. 

The form of the oscillation theorem used in this paper is that 
given by Kahn (1990). Different forms of this theorem were 
considered by mathematicians depending on specific differ- 
ential equations and the applied boundary conditions. The 
basic idea was originally introduced by Sturm (1836), who 
considered a second-order ordinary differential equation writ- 
ten in its standard form and established a comparison crite- 
rion, which allowed him to determine whether solutions to 
the equation were oscillatory or not, without formally solv- 
ing the equation. Another comparison criterion was devel- 
oped by Kneser (1893) and later generalized by Fite (1918), 
Hille (1948), Wintner (1949, 1957), and many others. The 
most important generalization of Kneser's criterion was done 
by Leighton (1950, 1962). In more recent work, some os- 
cillation theorems were established for linear (e.g., Li & Yeh 
1995, 1996) and nonlinear (e.g., Li 1998; Lee et al. 2005; 
Tyagi 2009) second-order differential equations. 



The oscillation theorem presented by Kahn ( 1 990), which is 
used in this paper, is equivalent to Kneser's criterion, which is 
also known as Kneser's theorem (e.g., Bohner & Unal 2005). 
First applications of this theorem to solar physics problems 
were done by Musielak & Moore (1995), who considered the 
propagation of linear Alfven waves in an isothermal solar at- 
mosphere. Then, Schmitz & Fleck (1998) used the theorem 
to establish criteria for acoustic wave propagation in the solar 
atmosphere. More recently the theorem was used by Schmitz 
& Fleck (2003), Musielak et al. (2006, 2007), Routh et al. 
(2007, 2010), and Hammer et al. (2010). 

In this paper, we use this method to derive the cutoff fre- 
quency for transverse waves propagating along a thin and 
non-isothermal magnetic flux tube embedded in the solar at- 
mosphere. The effects of temperature gradients on the cut- 
off frequency are studied for several power-law temperature 
models as well as for the reference mean solar atmosphere 
model C given by Vernazza, Avrett & Loeser (1981). The 
height dependence of the cutoff frequency in these models is 
calculated, and it is shown that the value of this cutoff at a 
given atmospheric height determines the frequency that trans- 
verse tube waves must exceed in order to be propagating at 
this height. The results are compared to those previously ob- 
tained for a thin and isothermal magnetic flux tube. We also 
briefly discuss implications of our results for the energy and 
momentum balance of the solar atmosphere. 

Our paper is organized as follows: the governing equations 
and derivation of our wave equations are given in Section 2; 
the global cutoff frequency for linear transverse waves prop- 
agating along a thin and isothermal magnetic flux tube is ob- 
tained in Section 3; the local cutoff frequency for the wave 
propagation along a thin and non-isothermal magnetic flux 
tube is derived in Section 4; the local cutoff frequencies for 
different power-law temperature distributions and for the VAL 
C solar atmosphere model are presented and discussed in Sec- 
tions 5 and 6, respectively; and conclusions are given in Sec- 
tion 7. 

2. WAVE EQUATIONS 

We consider a thin and non-isothermal magnetic flux tube 
that is embedded in the solar atmosphere. The tube axis is 
assumed to have a circular cross-section and an unperturbed 
tube axis oriented vertically along the z-axis, so that gravity 
g = —gt The tube density, pressure, temperature, and mag- 
netic field are respectively given by po = po(z), po = Po(z), 
To = Tq(z), and Bq = Bq(z). Moreover, the density, pressure and 
temperature of the external non-magnetic (B e = 0) atmosphere 
are represented by p e = p e (z), p e = p e (z) and T e = T e (z), respec- 
tively. We assume that the tube is in thermal equilibrium with 
its surroundings, which means that at each atmospheric height 
To(z) = T e (z). 

In the thin flux tube approximation (e.g., Ferriz-Mas et al. 
1989), the horizontal pressure balance for the tube becomes 

P0+^=Pe- (1) 

The equation for transverse linear oscillations of the tube 
was originally derived by Spruit (1981), who obtained 

9 2 t 2 . ,d 2 £ p -pedt 

-^2~C k (z) 77 -j-g ; -7- =0, (2) 

Ot dz 2 p0 + Pe OZ 

where £ is the horizontal displacement of the tube, and is 
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the characteristic wave speed given by 

Bo 



Ck(z) ; 



y/47T(po + p e ) 



(3) 



The relationship between the displacement £ and the corre- 
sponding magnetic field perturbation b x was obtained by Stix 
(1991) and Bogdan et al. (1996), and it can be written as 



dz 

Defining the velocity perturbation 

dt ' 

and using 



v,= ■ 



g 



PO + Pe 



2H 



(4) 



(5) 



(6) 



where H = c\/^g is the pressure scale height, with 7 being the 
ratio of specific heats and cs being the speed of sound given by 
c s = v7Po/po> we write Equation (0 in the following form 



d 2 v x 

dt 2 



-4(z) 



d 2 



cj(z) dv x 

dz 2 T 2H(z) dz 



(7) 



Then we combine Equations (0 and (0, use dBo/dz = 
-Bq/2H, and obtain 



d 2 b x 
dt 2 



-cliz) 



d 2 b x 
dz 2 



-cliz) 



( dc k (z) \ 1_ 

J 2H(z 



c k (z) V dz 



dz 



(8) 

The derived wave equations describe the propagation of linear 
transverse waves along a thin and non-isothermal magnetic 
flux tube. Since the wave equations for v x and b x are differ- 
ent, there is a phase difference between the wave variables. 
Physical consequences of the existence of this phase shift will 
be explored in the next sections. 

3. GLOBAL CUTOFF FREQUENCY 

We now consider the special case of a thin and isothermal 
magnetic flux tube by taking 7b = const, T e = const and 7b = T e , 
which gives c k = const and H = const. This allows us to write 
the wave equations (see Equations (0 and (0) as 



d 2 v x 2 9 2 v.x c\ dv x 
dt 2 ° k dz 2 2H dz 



and 



d 2 b x 2 d 2 b x 



dt 2 



dz 2 



2H dz 



(9) 



(10) 



Comparison of the above wave equations shows that they are 
different, which means that the behavior of the wave variables 
v x and b x is different even for a thin and isothermal flux tube. 

The wave equation for v x was originally obtained by Spruit 
(1981, 1982) and he showed that the propagation of trans- 
verse waves along a thin and isothermal magnetic flux tube is 
affected by a cutoff frequency. A new result presented in this 
paper is that we also derive the wave equation for b x , show that 
it is different than the one obtained for v x , and demonstrate 
how to obtain the cutoff frequency for both wave variables. 

To obtain the cutoff frequency, we follow Musielak & Ulm- 
schneider (2001) and cast Equations (0 and ( fTOb into their 
standard (or Klein-Gordon) forms (e.g., Roberts 1981; Rae 



& Roberts 1982). The required transformations are v x (z,t) = 

-1/4 1/4 

v\(z,t)p fe) and b x (z,i) = b\{z,t)p Q (z), and the resulting 
wave equations can be written as 



0» 
dt 2 



d 2 
k dz 2 



[vi(z,f),*i(z,0] =0, (11) 



where 



477 



(12) 



is called here Spruit's cutoff frequency. Since c k = const and 
77 = const in the isothermal case, J\ is also constant, which 
means that it is a global quantity. 

With the coefficients of Equation ( fTTT i being constant, we 
can make Fourier transforms in time and space, and derive 
the global dispersion relation: (u 2 - il£) = k 2 ^, where ui is 
the wave frequency and k = k z is the wave vector along the 
tube axis. Based on this dispersion relation, ilk is the global 
cutoff frequency for transverse tube waves. According to the 
dispersion relation, the waves are propagating when u> > il k 
and k is real, and non-propagating when either u> = f\ with 
k = or u> < fik with k being imaginary; in the latter case, the 
waves are called evanescent waves. 

According to Musielak et al. (2006, 2007) and Petukhov 
(2006), gradients of restoring forces, which can introduce 
characteristic scale heights, are responsible for the origin of 
cutoff frequencies for different waves. Using this argument, 
Musielak et al. (2006) argued about the origin of the acoustic 
cutoff frequency and Musielak et al. (2007) showed that the 
propagation of torsional tube waves is cutoff-free. Moreover, 
Petukhov (2006) demonstrated that horizontal, non-uniform 
but potential magnetic fields have no influence on a cutoff fre- 
quency for magneto-acoustic waves and, as a result, the cutoff 
is identical to Lamb's cutoff frequency, f^s = C&/2H, which 
was obtained by Lamb (1908, 1975) for acoustic waves prop- 
agating in a stratified and isothermal medium. 

The results obtained here for transverse tube waves clearly 
show that the vertical and non-uniform (exponentially diverg- 
ing) magnetic field of a thin and isothermal flux tube is re- 
sponsible for the origin of Spruit's cutoff frequency. Actually, 
the existence of this cutoff is caused by the horizontal pressure 
balance, which relates the tube magnetic pressure to the gas 
pressure. Since the latter introduces the pressure scale height 
H that also determines the rate with which the field diverges 
with the atmospheric height, the tube magnetic field and its 
characteristic scale height are the main physical reasons for 
the existence of Spruit's cutoff frequency. 

Finally, we may relate the Spruit cutoff ilk to the Lamb 
cutoff fis by writing 



ilk = — — f2s 

2c s 



(13) 



which shows that the ratio of Ck to cs determines how much 
ilk differs from lis- For the special case Ck = cs, we have 
flk = fls/2. Hence, flk can be much larger than fls when 
c k > > cs and much smaller when Ck < < cs . 

4. LOCAL CUTOFF FREQUENCY 

We now consider a thin but non-isothermal magnetic flux 
tube whose internal temperature distribution is To = 7b(z), and 
assume that at each atmospheric height the tube is in thermal 
equilibrium with its surroundings, which means that Tq(z) = 
T e (z). The presence of temperature gradients makes both c k 
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and H to be functions of z, and this leads to wave equations for 
v x and b K that have non-constant coefficients (see Equations 
© and dTOl). 

Since these wave equations cannot be solved by using 
Fourier transforms in space, a different method to obtain cut- 
off frequencies is required. Such a method was originally de- 
veloped by Musielak et al. (2006) for acoustic waves prop- 
agating in non-isothermal media. Routh et al. (2007, 2010) 
and Hammer et al. (2010) used the method to determine cut- 
off frequencies for the propagation of torsional waves in thick 
(isothermal) and thin (non-isothermal) magnetic flux tubes. 
Here, this method will be used to derive local cutoff frequen- 
cies for transverse waves propagating along a thin flux tube 
with different temperature profiles. 

The method is based on the oscillation theorem given 
by Kahn (1990), which follows the original work of Sturm 
(1836) and Kneser (1893). Actually, there are numerous oscil- 
lation theorems developed by mathematicians for various dif- 
ferential equations and their boundary conditions (e.g., Swan- 
son 1968; Teschl 2011, and references therein). However, 
most of these theorems cannot be directly applied to Euler's 
equation (e.g., Wong 1996; Aghajani & Roomi 2012), which 
is used in this paper. Fortunately, solutions to Euler's equa- 
tions are well-known (e.g., Murphy 1960), and they form the 
basis for our method used in this paper. 

4. 1 . Transformed wave equations 
The method begins with introducing a new variable 

dz 



dr ■ 



(14) 



The physical meaning of this variable becomes obvious after 
both sides of the above equation are integrated (see below). 
Then, r(z) is the actual wave travel time between a height at 
which a wave source is located and a given height z along the 
axis of the magnetic flux tube. 

We express Equations (|7]i and (0 in terms of the variable r 
and obtain the following transformed wave equations 



and 



"5 2 


d 2 




)-} 


dt 2 


dr 2 + 


W 2#, 


1 dr 


' d 2 


•-"( 




d ' 


dt 2 


dr 2 \ 


,c k 2HJ 


dr. 



v x (T,t)=Q (15) 



b x (T,t) = 0, 

where c k = Ck(r), c' k = dc k /dr and H = H(j). 



(16) 



4.2. Standard wave equations 

To convert the transformed wave equations into their stan- 
dard forms, we use 



v x (r,t) = v(T,t)exp 

and 

b x (T,t) = b(T 7 t)exp 
and obtain 



1 

+ 2 



et 2H 



~~ 1 di~ 



4 + _Ck 

c k 2H 



rr: df 



dr 2 



d 

~Q^2 ~ + ^ci-,y( T ) 



v(r,f) = 



(17) 



(18) 



(19) 



and 
where 

and 

a 2 ,., 



dt 2 ' 



dr 2 



Kr,t) = 0, 



n 2 (T) = n 2 (T)+-- 



1 r" 

2 c k 



k+ c k H> 



AH H 



(20) 



(21) 



,(r): 



14' 
+ - — - 



2 c k 



^ + ^, (22) 
AH H 2H ' 



with c(.' 



d 2 Ck/dr 2 and H' = dH /dr. The frequencies £l CI V 
and £l cr ^ are known as the critical frequencies (Musielak et 
al. 1992, 2006). 



4.3. Turning-point frequencies 
We make the Fourier transform in time [v(T,t),b(r,t)] = 
[v(T),b(T)]e~' u ' , where uj is the wave frequency, and write 
Equations ( fT9l and d20b as 

d 2 



dr 2 



■+UJ 



and 



y(r) = 



b(r) = 



(23) 



(24) 



Using the oscillation theorem given by Kahn (1990) and 
comparing the above equations to Euler's equation (e.g., Mur- 
phy 1960), we obtain the following conditions for the wave 
propagation 

1 



and 



nl v (T)> 



At 2 



■n 



i 



(25) 



(26) 



Note that applications of the oscillation theorem to other wave 
propagation problems relevant to solar physics were previ- 
ously considered by Musielak & Moore (1995) and Schmitz 
& Fleck (1998); see also Routh et al. (2010) for more re- 
cent results. The idea of using Euler's equation to determine 
turning-point frequencies was first introduced by Musielak & 
Moore (1995). 

We follow Musielak et al. (2006) and define the turning- 
point frequencies as 



O 



1 



k v(T ) = f£ iV ( T ) + _ 



and 



(27) 



(28) 



where r is the actual wave propagation time (see Routh et al. 
2010) and is given by 



T(Z) = 



dz 
c k (z) 



(29) 



with tq being an integration constant to be evaluated when 
flux tube models are specified (see Sections 6 and 7). 

It is important to point out that the turning-point frequen- 
cies were obtained by making the derived wave equations (see 
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Equations d23l l and d24i i) equivalent to Euler's equation, so 
that the well-known solutions of the latter can be directly ap- 
plied to these wave equations. 

4.4. Converting r into z 

Having obtained the turning-point frequencies (see Equa- 
tions d27| i and (l28li), we now express them in terms of the z 
variable, which requires using Equation (f29) to convert r to z. 

We begin with the critical frequencies £l a ,v(T~) and fi C r,i(T) 
given by Equations (l2TT i and d22i i. Using the following ex- 
pressions 

1 dc k dc k 
Ck dr dz ' 



(30) 



1 d z c k 
Ck dr 2 



dz 1 



dc k 
dz 



and 



we obtain 



]_dH 

Ck dr 



dH 

dz 



(31) 



(32) 



1+4| — 

dz 



n 2 a .jz) = n 2 k (z) 

and 



1 d Ck 1 / dc k 
'2 CiL ~dz T + 4 \~dz~ 



1-4 



dH 
~dz 



Ck dc k 
2H dz 



1 d 2 c k 1 / afck 
+ 2 Ck "5?" + 4 I dz" 



(33) 



(34) 



We now use Equation ((29) to express the conditions for 
wave propagation and the turning-point frequencies as func- 
tions of z, and obtain 



and 



[uj 2 -n 2 cl . tV (z)]> 1 - 



[w 2 -Q^,(z)] > - 



' dz 
z dz 



1-2 



+ T C 



The turning-point frequencies are 



and 



n( p jz) = nUz)+- 



nl b (z) = nl/z)+- 



z dz 

z dz 
Ck(z) 



1-2 



(35) 



(36) 



(37) 



(38) 



We shall use the turning-point frequencies given above to 
determine the cutoff frequency for transverse tube waves. 

4.5. Cutoff frequency 

We follow Musielak et al. (2006) and take the larger of the 
two turning-point frequencies to be the cutoff frequency Q cut , 

n cut (z) = max[0^ )V (z), Qtp,b(z)] ■ (39) 

Our selection of fi cut is physically justified by the fact 
that in order to have propagating transverse tube waves, the 



10' 
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Figure 1. Temperature is plotted versus the distance ratio z/zo for the power- 
law temperature models with m ranging from 1 to 5. 



wave frequency us must always be higher than any turning- 
point frequency. In other words, the choice guarantees that 
propagating wave solutions are obtained for both wave vari- 
ables, and that the cutoff frequency does separate the prop- 
agating and non-propagating wave solutions (see also Ham- 
mer et al. 2010). Hence, the condition for propagating waves 
is u > ficut- Based on our definition of the turning-point 
frequencies, the condition for non-propagating (evanescent) 
waves is to < £! cut . 

The above results show that the cutoff frequency can only 
be determined when a flux tube model is specified, so that 
the turning-point frequencies can be obtained (see Sections 5 
and 6). In addition, one must keep in mind that the condi- 
tions given by Equation d39l must be checked at each height 
because in some regions along the tube il lp v could be larger 
than Sltp./j, and the opposite could be true in other regions (see 
Section 6). 

5. MODELS WITH POWER-LAW TEMPERATURE 
DISTRIBUTIONS 

Let us consider the following temperature distribution in- 
side the tube 



T (z) = T {m C 



(40) 



where £ = z/zo is the distance ratio, with zo being a fixed 
height in the model, Too is the temperature at zo, and m can 
be any real number. Note that in all power-law models dis- 
cussed below a wave source is assumed to be located at £ = 1, 
which means that in all calculations £ > 1. In addition, for all 
models we take zo = 10 km, Too = 5000 K, Cko = 10 km/s, and 
for gravity we take its solar value. The resulting temperature 
distributions for m being a positive integer that ranges from 1 
to 5 are presented in Figure 1 . 

5.1. Case of m = 1 

To describe the process of deriving a local cutoff frequency, 
we begin with the simplest case of m = 1, which corresponds 
to the temperature varying linearly with £. We calculate po, 
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Figure 2. The normalized cutoff frequency S7 cut / f!k0 is plotted versus the 
distance ratio z/zo for the power-law temperature models with m ranging 
from 1 to 5. 



Pq, Bq and Ck as functions of £, and use Equation (|29l to obtain 

T(0 = 2—e l2 + rc, (41) 

where Cko is the value of Ck at zo and rc is the integration 
constant. To determine this constant, we assume that t(£ = 
1) = r = zo/cko, which gives r c = -r . 

We express the turning-point frequencies given by Equa- 
tions d37l > and d38l in terms of £. This gives 



kO 



(42) 



where Oko = Cko/4#oo, with Z/oo being the scale height at zo, 
and 

gl(0= 7^7T# 7^ (43) 



and 



k() 



(2^V2-1)2 



zo 



(44) 



The difference between the turning-point frequencies is 



^,(0-^(0 



"kO 



= 4 



#00 

zo 



ZO 



(45) 



which shows that 17 2 p v (£) is larger than f2 2 p Hence, O t p,i 
becomes the local cutoff frequency, so we write ri C ut(0 = 
^tp,v(0 and 



fl cut (0 = a 



k() 



1+4^ + 

zo 



#00 

ZO 



(3 + *i(0) 



1/2 



-1/2 



(46) 

The local cutoff frequency O cut is plotted as a function of £ in 
Figure 2. It is seen that the cutoff frequency decreases with 
the atmospheric height in the model with m=l. 



5.2. Case of m = 2 
In this case, t(£) is given by 

t(0 = — Int + Tc, 

CkO 



(47) 



where tq is the integration constant determined from the as- 
sumption that t(£ = 1) = ro = zo/cko- This gives r c = to. The 
turning-point frequencies are 



r 2 + 8 ^ r . +4 te) (l +ft(G ) 



Z() 



Z() 



and 



where 



(^) 



*2(0 = 



1 



(1 + ZnO 2 



(48) 



(49) 



(50) 



Here f2tp,v(0 is always larger than il tp .b(0- Thus, we choose 
f2 tp iV as the local cutoff frequency and write n cut (0 = Otp jV (0, 
or 



^cut(0 = ^k 



r 2 +8 ^22 r i +4 



zo 



ZO 



H +82(0) 



1/2 



(51) 

This cutoff frequency is plotted as a function of £ in Figure 
2, which shows that the cutoff remains almost constant in the 
model with m = 2. 

5.3. Cases with m > 2 
In this general case of m > 2, we obtain 



, t , zo 
t(0 = — 



c M \2-m 



(52) 



where the integration constant tq is evaluated by taking t(£ = 
1) = To = Zo/cko; note that our choice of to gives the same 
physical parameters at z = Zo for all the power-law models. 
After evaluating tc, we calculate 



where 



+ „„4-,„)('^Vr- ! 

m 1 V zo / 



(53) 



(54) 



and 



^(0 = ^„ 



»J / V zo / 



(55) 
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We now calculate the difference between these turning- 
point frequencies and obtain 

"U>-nU> =4m(m . 2) f jM V-^gy . 

O k0 V zo J 20 

(56) 

Since in this model Hqq > Zo, m > 2, and £ > 1, fi t 2 p 6 is larger 
than f2^, v . This indicates that the local cutoff frequency is 

^cut(0 = ^tp,fc(0, or 



k() 



r m +'«(3m-4) 



#oo 

20 



e 



■m-1 



ml \ zo J 



1/2 



(57) 



The cutoff frequency calculated for the power-law tempera- 
ture models with m = 3, 4 and 5 is plotted versus the dis- 
tance ratio in Figure 2. It is seen that this cutoff frequency 
always increases with the atmospheric height in the models 
with m > 2 and that its increase is much faster for higher val- 
ues of m. 

5.4. Discussion 

The effects of different temperature gradients on the cutoff 
frequency for transverse tube waves are presented in Figure 
2. Since the cutoff frequency is a local quantity, its value at 
a given atmospheric height determines the frequency that the 
waves must exceed in order to be propagating waves at this 
height. Our results demonstrate that the conditions for wave 
propagation strongly depend on the temperature profiles. 

If the temperature increases linearly with height ( m =1), the 
cutoff frequency starts with its maximum at z = Zo and then 
decreases with height. For the temperature model with m = 2, 
the cutoff frequency remains practically constant with height. 
However, the cutoff frequency always increases with height 
in the temperature models with m > 3; the higher the value of 
m, the steeper an increase of the local cutoff frequency with 
height is observed. 

The main purpose of using the power-law temperature mod- 
els was to demonstrate the dependence of the local cutoff fre- 
quency on the increasing steepness of the temperature models. 
Obviously, the power-law models do not properly describe the 
temperature gradients in the solar atmosphere. Therefore, we 
now consider a more realistic model of the solar atmosphere. 

6. VAL MODEL OF THE SOLAR ATMOSPHERE 

We assume that a thin flux tube is embedded in the refer- 
ence mean solar atmosphere model VAL-C (Vernazza et al. 
1981). In this model, the height z = km corresponds to 
unity optical depth, where the temperature is 6420 K. Beyond 
the temperature minimum r m! „ = 4170 K, which is located at 
Z = 5 15 km, the model extends through the chromosphere well 
into the transition region to the corona, up to a temperature of 
4.47 x 10 5 K. At the base of the transition region the VAL-C 
model exhibits a plateau with a temperature around 2 x 10 4 
K, caused by Lya emission. This plateau was later shown to 
disappear if the physics of the transition region is treated more 
realistically (Fontenla et al. 1990). Therefore we restrict our- 
selves to the height range from z = Zo = km (optical depth 
unity) to the just below the plateau (z = 2 1 13 km). 



To calculate the characteristic wave speed Ck as a function 
of z, we have to know the magnetic field Bo(z), which is not 
given in the VAL-C model. Since the tube is assumed to be 
in in thermal equilibrium with its surroundings at each atmo- 
spheric height, we use the local value of the pressure scale 
height, which is the same outside and inside the tube, to eval- 
uate Bq(z); the calculation begins with an assumed surface 
magnetic field Bq(zq) = 1500 G (see Equation (Q])). Having ob- 
tained the distribution of the tube magnetic field with height, 
we then calculate the gas pressure and density inside the tube, 
and the characteristic wave speed Ck(z) (see Equation (O). In 
the lower panel of Figure 3 we plot Ck versus the atmospheric 
height in the model; for comparison we also plot the sound 
speed cs- The results show that Ck is smaller than cs in almost 
the entire model, except in the upper chromosphere and lower 
transition region, where Ck becomes comparable to, or even 
slightly larger than, cs. 

In order to calculate the wave propagation time t, we use 
Equation d29l ). For the integration constant tc = t(z = 0) = to 
we choose a value of 100 s, thus assuming that the waves have 
traveled for 100 s before they reach the base zo of the model. 
For typical wave speeds (cf. lower panel of Figure 3) this 
means that the waves are generated about two scale heights 
below zo- 

To determine the cutoff frequency Sl C ut(z) according to Sec- 
tion 4.5, we calculate both turning-point frequencies Qtp,v(z) 
and ilip.b(z) and select the larger one as f2 cu t(z) (see Equations 
d33l . ( 134b . and d37]i-(l39l). To perform these calculations, we 
must evaluate the first derivative of Ck and H, and the second 
derivative of Ck. These calculations have to be done numeri- 
cally because the model consists of tabulated data. The upper 
panel of Figure 3 shows the resulting cutoff frequency f2 C ut(z) 
as well as Spruit's cutoff frequency £lk(z) = c^(z)/4-H(z) (cf. 
Equation (fT2l) ). where the latter is again treated here as a lo- 
cal, height-dependent quantity. 

A striking property of Figure 3 is that the cutoff frequency 
for the nonisothermal case, fi C ut(z), is wiggly. This is mostly 
caused by the dependence of the turning-point frequencies on 
the second derivative of the tube speed Ck (cf. Equations d33l 
-d34ll). Since the VAL-C model is specified by a number of 
tabulated distinct points, the curves in Figure 3 must be plot- 
ted by finding interpolating or fitting functions through these 
points and determining the first derivative of the scale height, 
and the first and second derivatives of the tube speed, in order 
to calculate the terms in Equations (|33T>-(|341>. Such deriva- 
tives of interpolated curves tend to show a wiggly behavior. 
As shown in the middle and lower panels of Figure 3, these 
wiggles are stronger near the temperature minimum and in 
particular at the foot of the transition region, where the gradi- 
ents change rapidly (see also Fawzy & Musielak 2012). 

One could avoid the wiggles by approximating the data by 
simple functions with known smooth derivatives. There are 
also automatic algorithms that can smooth out such wiggles 
more or less completely (e.g., Chartrand 201 1). However, we 
show them here explicitly in order to illustrate the dependence 
of the cutoff frequency, via these derivatives, on the nature of 
the solar atmosphere, which in reality is less homogeneous 
than the VAL-C model and in addition temporally variable. 
Therefore, we chose the well-proven algorithm by Reinsch 
(1967), which fits the data points by a natural cubic spline, 
thus giving consistent first and second derivatives. A mod- 
erate amount of smoothing can be applied, so that the spline 
needs not go exactly through the points. 

Since the value of f2 cut is different at each atmospheric 



8 



Routh, Musielak, & Hammer 



0.2 



0.15 — 



0.1 — 



* 0.05 — 



U 



X 



— 



cut 
T 




18 

16 

14 

12 

10 

8 

6 

4 



500 1000 1500 2000 

Height z [km] 




500 1000 1500 

Height z [km] 



2000 



„ 20 



15 



10 



5 — 



l — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — r 



dc,,/dz 




j I I I I I I I I I I I I I I I I I I I I i 



0.4 



0.3 



0.2 



0.1 



— 



500 



1000 
Height z [km] 



1500 



2000 



-0.1 



Figure 3. The upper panel shows the height dependence of the cutoff frequency f2 cut for the VAL-C model, for which the temperature stratification is shown as 
a dotted line scaled with the right y axis. For comparison, Spruit's isothermal cutoff frequency = c\/4H is plotted by using local values of the tube speed cjj 
and the scale height H. The middle and lower panels, respectively, show H and c\ as well as the sound speed c's, with scales on the left y axes. The values of H 
and Ck at the tabulated points of the VAL-C model are marked with circles. The derivatives of the spline curves through these points are shown with scales on the 
right y axes. 



height, the frequency w of a transverse tube wave must be 
higher than the cutoff at a given height z in order to be prop- 
agating at this height. Overall, the new cutoff frequency il cut 
is seen to be slightly higher than ilk throughout the chromo- 
sphere. This is mostly caused by the fact that the chromo- 
spheric temperature increases with height, leading to a posi- 
tive first derivative of the tube speed c^, which increases both 
turning-point frequencies (cf. Equations (l3Jt and (|34| |) and 
thus also their maximum, the cutoff frequency i7 cut (Equation 
03). 



Another difference between the Spruit cutoff and our new 
result lies in the behavior at the foot of the transition region. 
Here ilk decreases with height (which becomes more evident 
if one continues the calculation to higher levels of the VAL- 
C model), while il cut increases rapidly as we approach the 
steep temperature rise in the transition region. Formally this 
means that a wave with arbitrary frequency tends to lie above 
ilk at the transition region (implying propagation in the case 
of an isothermal treatment of the cutoff), but below our value 
(implying reflection off the temperature gradient in our non- 
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isothermal treatment). This may have implications for the 
contribution of such transverse waves to the heating and dy- 
namics of the overlying corona. However, we put not much 
weight on such differences in the high chromosphere and tran- 
sition region, fow two reasons. First, we have not provided a 
complete study of the reflection properties of the transition 
region, which would require a more complex analysis of the 
complete solar atmosphere. And second, a basic approxima- 
tion underlying all these theories breaks down much earlier, 
already in the mid chromosphere - namely, that the flux tubes 
can be approximated as thin. The thin flux tube approxima- 
tion, which is the basis for the results presented in this paper 
and for Spruit's results as well, is no longer valid at those 
heights. Thus the results obtained for the upper atmosphere 
must be taken with caution. We have extended the calcula- 
tions up to the foot of the transition region in order to deter- 
mine the approximate behavior of the cutoffs in these layers, 
even though this is beyond the formal limit of validity of our 
analytical results, because the presented results demonstrate 
that the temperature gradient in the upper solar chromosphere 
and in the solar transition region will have a major influence 
on the propagation of transverse tube waves. 

The increase of temperature with height requires atmo- 
spheric heating, which is typically identified with acoustic- 
gravity and flux tube waves, including transverse tube waves, 
or with phenomena related to magnetic reconnection (e.g., 
Priest 1982; Narain & Ulmschneider 1996; Ulmschneider & 
Musielak 2003). Our results presented in Figure 3 give con- 
straints on the range of frequencies of transverse tube waves 
that are propagating in different parts of the solar atmosphere. 

Our results show that in the upper photosphere and lower 

chromosphere f^ t ot ft O.Olrad/s, or the cutoff period is 500 
s, and that in the middle and upper chromosphere 0, c J^ om ft 
0.02rad/s, which corresponds to the cutoff period of 350 
s. As already mentioned in Section 1, Fujimura & Tsuneta 
(2009) observed transverse tube waves with periods of 240- 
540 s, or the corresponding frequency range u ft (0.01 - 
0.03)rad/s, in intergranular magnetic elements. Comparison 
of these observational results to our theoretical cutoff frequen- 
cies shows that most of the observed waves are freely propa- 
gating in the solar photosphere and lower chromosphere be- 
cause id > ricut" 1 , however, in the middle and upper parts of the 
solar chromosphere the waves with frequencies us < r^™™ 
are not propagating. 

Now, the high frequency waves with periods of 45 s, or 
lo ft 0.1rad/s, observed by Okamoto & De Pontieu (2011) 
along spicules do satisfy the condition uj > fi* om and they 
are obviously freely propagating waves in the solar chromo- 
sphere. Similarly, the high frequency oscillations with peri- 
ods ranging from 30 s to 180 s, which corresponds to the fre- 
quency range ui ft (0.03- 0.2) rad/s, detected by Yurchyshyn 
et al. (2012) in type II spicules also satisfy the condition 
uj > ^cm™' an d therefore they are freely propagating waves 
in the solar chromosphere. 

The above constraints can be used to assess the role of 
transverse tube waves in the heating and dynamics of the so- 
lar atmosphere. As discussed by Hasan & Kalkofen (1999), 
Musielak & Ulmschneider (2003), and Hasan (2003), trans- 
verse tube waves may be responsible for the excitation of so- 
lar atmospheric oscillations observed in magnetically active 
regions near sunspots. The results obtained in this paper can 
be used to determine the natural frequency of the solar atmo- 
sphere inside thin and non-isothermal magnetic flux tubes and 



the effects of temperature gradients on solar atmospheric os- 
cillations. 



7. CONCLUSIONS 

Our method to determine the cutoff frequency for trans- 
verse waves propagating along a thin and non-isothermal flux 
tube requires that integral transformations are used to cast the 
wave equations for both wave variables in the standard forms. 
Then the conditions for the existence of propagating wave so- 
lutions are established by using the oscillation theorem. The 
theorem is also used to obtain the turning-point frequency for 
each wave variable. The larger among the two turning-point 
frequencies is selected as the cutoff frequency. 

To study the effects of temperature gradients on the cutoff 
frequency, we used different power-law temperature distribu- 
tions. For the temperature that increases linearly with height, 
the cutoff frequency is large at the bottom of the model and 
then decreases with height. For the temperature profile with 
the power of 2, the cutoff frequency remains practically con- 
stant with height. However, for powers higher than 2, the cut- 
off frequency always increases with height; and the higher the 
power, the steeper the increase of the cutoff frequency with 
height. 

We also calculated the cutoff frequency fl cut as a func- 
tion of height z in the VAL-C model and compared it to 
Spruit's global cutoff frequency fis that was treated as a 
height-dependent quantity. The comparison shows that f2 cut 
exceeds £1$, and that the differences are especially prominent 
in the upper parts of the model, where however the thin flux 
tube approximation is not valid any longer. The differences in 
the lower atmosphere, where the thin flux tube approximation 
is valid, may be important for the energy carried by transverse 
tube waves from the solar convection zone, where the waves 
are generated, to the overlying solar atmosphere, where the 
wave energy is deposited. 

The cutoff frequency O cut calculated for the VAL-C model 
gives constraints on the range of frequencies of transverse 
tube waves that are propagating in different parts of the so- 
lar atmosphere. We compared f2 cut to the range of frequencies 
determined from observations by Fujimura & Tsuneta (2009), 
Okamoto & De Pontieu (201 1), and Yurchyshyn et al. (2012). 
Based on this comparison, we established that most of the 
waves observed by Fujimura & Tsuneta (2009) in intergranu- 
lar magnetic elements are freely propagating in the solar pho- 
tosphere and lower chromosphere, however, in the middle and 
upper solar chromosphere the lower frequency part of these 
waves corresponds to non-propagating waves. On the other 
hand, all the waves detected by Okamoto & De Pontieu (201 1) 
and Yurchyshyn et al. (2012) in spicules are freely propagat- 
ing waves in the solar chromosphere. 

These constraints are important for understanding the role 
played by transverse tube waves in the heating of the solar at- 
mosphere and in the acceleration of plasma, e.g. in spicules 
and the solar wind. Our results can be used to determine the 
natural frequency of the solar atmosphere inside thin and non- 
isothermal flux tubes and to study the effects of the tempera- 
ture stratification on the excitation of atmospheric oscillations 
inside solar magnetic flux tubes. 
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